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Abstract 

We set up a formalism that can be used to calculate the power spectrum of the 
curvature perturbations produced during inflation up to arbitrary order in the slow- 
roll expansion, and explicitly calculate the power spectrum and spectral index up to 
second-order corrections. 
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1 Introduction 



Curvature perturbations produced during inflation are thought to be the origin of the inho- 
mogeneities necessary for galaxy formation and other large-scale structure. It will soon be 
possible to measure the power spectrum of these inhomogeneities very accurately by means 
of the MAP [0] and Planck @] satellites and the SDSS program ||. For example, Planck 
is supposed to be able to measure the power spectrum of the cosmic microwave background 
radiation with an error of a few percent ||, and the spectral index with an error of less than 
one percent [|J . Thus, it is very important to calculate the power spectrum precisely in order 
to get the maximum benefit from these observations. 

The power spectrum and spectral index of the curvature perturbations produced during 
inflation have been calculated up to next-to-leading order, or first-order corrections, in the 
slow-roll approximation H ||. Given the current observational constraints on the spectral 
index, these results for the power spectrum and spectral index are expected to have errors 
of order a few percent or less and a few tenths of a percent or less, respectively. These errors 
are comparable with the expected errors in the planned observations discussed above. It is 
thus desirable to improve the accuracy of these calculations to ensure that the calculational 
errors are small compared with the observational errors. See Refs. |7j for further discussion 
of the accuracy of, and previous efforts to go beyond, these results. 

In this paper, we set up a formalism that can be used to calculate the power spectrum of 
the curvature perturbations produced during inflation up to arbitrary order in the slow-roll 
expansion, and explicitly calculate the power spectrum and spectral index up to second- 
order corrections. To check our result, we compare it with the two known exact solutions 
||, power-law inflation and inflation near a maximum. Our results are only valid for a 
single component inflaton, but will be generalized to a multi-component inflaton H, ^] in a 
forthcoming paper. 



2 The Calculation 

2.1 Preliminaries 

The scalar linear perturbation to the homogeneous, isotropic background metric is generally 
expressed as || 

ds 2 = a 2 (r]) {(1 + 2A)drf L - 2^5 dx*drj - [(1 + 211)5^ + 2^8 jH T ] dx l dx j ) . (1) 

It is convenient to express the density perturbation in terms of the intrinsic curvature per- 
turbation of comoving hypersurfaces 1Z C [[UJ, given by 



K C = K- ^50 (2) 



during inflation, where 5(f) is the perturbation in the inflaton field. The power spectrum 
Vn c (k) is given by 

(^ c (x,r / ),^ c (y,r / ))= / — ^ '_ yu V^. (3) 
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The standard result for the power spectrum of the curvature perturbation to leading order 
in the slow-roll approximation is |Tl| 



f H\ 2 fH x 2 



and to first-order corrections is || 

VnM = [1 + (4a - 2)e x + 2a8 x \ 

where 



(4) 



aH=k 



H\ 2 fH x ~ 
2^ 



(5) 



aH=k 



El 



H 1 
IP ~ 2 \ H 



Si 



H<f> 



a = 2-ln2-7~ 0.729637 

and 7 is the Euler-Mascheroni constant, 7 ~ 0.577216. 

Our starting point is the effective action during inflation, 



(6) 
(7) 



S 



- l -R+ l -{d^f~V^) 



Defining 



: = — and ip = a I 5<fi — —TZ 
H \ H 



d 4 x. 



-ZRr 



(8) 



(9) 



the action for scalar perturbations is |T 

/ 1 / < i ~ \ 

S 



dip 
dr] 



dr] dx 3 



(10) 



where the conformal time rj is given by dr] = dt/a. Using the Fourier transform of ip, the 
equation of motion becomes 



d 2 f k + A 2 _ = () 



and <y?fc satisfies 



_J_ p ~ikri 
V2k 

A^z as — krj — > 0. 



( ""' as •/.•// • dc 



(12) 



The first asymptotic solution corresponds to the usual flat space vacuum on scales much 
smaller than the Hubble distance, and the second is the growing mode on scales much larger 
than the Hubble distance. The power spectrum is given by 





lim 


<-Pk 


2 k 3 
= \A k \ 


1 2vr 2 J 


— krj— >0 


Z 


2vr 2l k \ 



(13) 
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2.2 The Formalism of the Calculation 

Our task is to solve E g. flTT|) with the boundary conditions Eq. ([12]) to eventually calculate 
Ak- We define y = v2kfk and x = —ki]. Then, the equation of motion for y and the 
asymptotic solutions are 

d 2 y ( ld 2 z\ . A , 



(15) 



dx 2 \ z dx 2 / 
and 

r °} x as x — > oo 



2fc AfcZ as a; — >• 0. 
Now, we can choose the ansatz that z takes the form 



x 

Then we have 

1 d 2 z 2 I 



z = -fQnx). (16) 



z dx 2 x 2 x 2 



+ ^gQnx), (17) 



where 



-3/' + /" 

and the equation of motion is 

d 2 y ( 2 \ 1 



;is) 



/ > i 1 — 2 h = ^^( ma; )y- ( 19 ) 

dx z \ x z J x A 
The homogeneous solution with the correct asymptotic behaviour at x — > oo is 

y (x) = (l + e^. (20) 

Therefore, using Green's function, Eq. flTSP with the boundary condition Eq. ( |15|) can be 
written as the integral equation 

if 00 l 

y(x)=yo{x) + - du — g{\nu) y{u) [y* {u)y {x) - yt(x)y (u)] . (21) 

^ J X u 

2.3 The Slow-Roll Expansion 

We will solve Eq. (|21|) perturbatively. Our small parameters will be the slow-roll [|T3| , |5|] 
parameters ej and 5i of Eq. @ and 

A - 1 ^ f99^ 
^ = (22) 
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evaluated at horizon crossing {aH = k), which we assume to satisfy e\ < £ and \8 n \ < £ n 
for some small perturbation parameter £. Note that the subscript denotes the order in the 
slow-roll expansion. 
We can expand 



xz = /(In re) = ^j(hi 



x) 



(23) 



n=0 



where, as we will see, /„/ fo is of order n in the slow-roll expansion. This expansion is useful 
for exp(-l/f) < x < exp(l/£). 
Using 



x 



-krj — —k [ — ~ 4? (! + e i + 3e i + 2e i 5 i) 
y a aH 



(24) 



we can express the /„'s in terms of €\ and the 5 n evaluated at aH = k. To leading order in 
the slow-roll expansion 



Similarly, we can obtain 



d 2 xz 



(d\nx) 



x=l 



M (84 + 9e 1 S 1 + S 2 ) 



(25) 



A;0 
IP 

k$ 



(-2e 1 -5 1 -6e?-4e 1 5 1 ) 



/o * ^ (1 + ei + 5e? + 36^) 



(1 - ei - 4e? - 36x5!) 



-2 Cl - ^ - 4e? - 3ei<5 



so we get, up to the appropriate orders for this paper 

1 H 2 
fo 

k 

fo 

k 

fo 

Eqs. (PI) and (G3) give 



(26) 
(27) 



aH=k 
1 1 ai?=fc 



ie? + 9ei5i + (J 2 | a 



(28) 
(29) 
(30) 



(31) 



z — ' n! 

n=0 

where # n is of order n in the slow-roll expansion and, up to the appropriate orders, 

fi fa 

5i - -3^r + ^ (32) 



92 



fi\ _ ^k 

fo/ fo 
4 



(33) 



We can expand y as 



(34) 



ra=0 



where yo{x) is the homogeneous solution, Eq. (p0|) , and y n (x) is of order n in the slow-roll 
expansion. 

We can now solve Eq. (|2Tj) perturbatively by substituting Eqs. (|3T|) and ([34]) and equating 
terms of the same order. 



2.4 The Calculation up to Second-Order Corrections 



Substituting g = g\ and y = y into the right hand side of Eq. ([21~D, we can easily find the 
first order corrections to y as 



:9i 



du -^y (u) [yo(u)y (x) - y*(x)y {u)} 



ir 



;9i 



-2ix 1 e lx + y* Q (x) \ du u 1 e 



-l jiu 



(35) 



Now considering up to second-order corrections, substituting g = g\ + g 2 In x and y = 
y + yi into Eq. ( plf ) gives two terms of the same order, due to the product of g\ with y 1 and 
g 2 with y . Let us call them 2/21(2) and y 22 {x), respectively. Explicitly, 



2/21 [x ) 



and 



\9i 



1 2 
-g5l 



du ^yi{u) [yo(u)y {x) - ya{x)y {u)) 



'-x~ l e ix + 



+ iyo{x) 



--x- 1 + -)< — 
3 3 

du u- l e- 2iu 



du u~ x e 2iu 
dv v^e 2iv 



(36) 



2/22 (X ) 



:92 



:92 



du — \nu y (u) [yl(u)y Q (x) - yo(x)y (u)} 



u- 



x 



1/ 



+ 21ni y^x) / duvT^e 



-l„2iu 



+ %o(x) 



(i-u -u 1 In u e 



2iu 



(37) 



Adding together Eqs. (^), (|3"5D, (|3T)|) and fl3"T|) gives the perturbative solution of Eq. 
for small £|lnx| up to second-order corrections. It should be accurate for exp(l/£) ^> x ^> 
exp(— 1/£). Taking the limit x — > with ^ln(l/x) fixed and small, we obtain, after some 
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calculation, the asymptotic form for y up to second-order corrections 

, 2 



y{x) 



91 
3 

92 

+ — 
6 



111 



9l 
18 



a 



-a 
3 



7T 



a H — a V itt [ a -\ — 

3 12 V 3 



4 + — + Z7T I a - - 



»7 _^ 2 r 



i 

3 



7T 



+ ^ s 1 In x 



+ \\ g i-92\x-\\nxf. 



(3* 



The exact asymptotic form for y in the limit x — > is given by Eq. flTSp. Expanding this 
perturbatively as in Eq. ([2rJ) for small £ ln(l/x), i.e. for x in the range 1 ^> x ^> exp(— l/£), 
gives the asymptotic form for y up to second-order corrections 

y(x) -> v^^/o + V^fcA fc / a z^lnz + iv^LWa aT^mx) 2 . (39) 

Comparing this with Eq. d38|) , the coefficient of x _1 is the desired result because it will 
give Ak up to second-order corrections. The coefficients of a; -1 In a; and x~ 1 (lnx) 2 simply 
give the consistent asymptotic behaviour, that is, proportional to z. Therefore, substituting 
Eqs. ([$2|) and ( p3|) into Eq. (|38|), matching the coefficient of x~ x with that in Eq.(^), and 
using Eq. fljjD , the power spectrum is 



(2tt) 2 / 2 



h ( 2 5tt 2 
2a^ + 3a 2 -4 + — 
Jo V Vl 



Therefore, substituting Eqs 
H 4 



(2tt) 2 2 



1), (H) and © 
1 + (4a - 2)ei + 2a5i + (4a 2 



23 



3a 2 + 2a - 22 + 



29tt 5 



3a 2 - 4 + 



+ 



77^ 
3 

5tt 2 



-a z + 



7T 



12; /o 



(40) 



(41) 



7T 



12 y \ 12 

where the right hand side should be evaluated at ail = k. The numerical values of the 
coefficients in this equation are 0.918549, 1.459274, 2.158558, 4.907929, 1.709446, 0.290097, 
respectively. Note the curious cancellations ensuring that the coefficients are of order one. 
The spectral index 

dhxP Uc 



nnAk) = 1 + 



dink 



can be easily calculated, and we obtain the result 
n nc {k) = 1 - 4ei - 25i + (8a - 8)e 2 + (10a - 6)ei$i - 1ab\ + 2aS 2 



(42) 
(43) 



+ 
+ 



-16a 2 + 40a - 108 + 



28tt^ 



-31a 2 + 60a - 172 + 



199vr 
12 



e 2 A 



-3a 2 + 4a - 30 + 



13tt 2 



+ -7a 2 + 8a - 22 + 



31tt' 
12 



6l<5 2 + 



ei5 2 



-2a 2 + 8 



5vr 5 



si 



3a 2 - 8 + 



3tt ; 



7T 



6^2 4 ( -a z + — ) 5 3 



where again the right hand side should be evaluated at aH = k. The numerical values of 
the coefficients in this equation are -2.162903, 1.296372, -1.459274, 1.459274, 4.783868, 
18.945687, 3.397652, -1.289411, 5.606983, 0.999314, 0.290097, respectively. 



2.5 Checking against the Exact Solutions 



To check this result, we use the two exact solutions, power-law inflation |Tj] and inflation 
near a maximum |5j . Note that these two solutions can be derived from the Hankel function 
solution H which corresponds to all orders in g\ with g n = for n > 1. 

2.5.1 Power-law inflation 

In power-law inflation, the potential is 



V(<j>) = V exp 



From this, 



1 r 

ei = - , <Ji 
V 



1 

V 



5 2 = 25\ , S 3 = 6Sf 



and the power spectrum and the spectral index are 

pH 2 



H 4 



(2n) 2 (f) 2 



P 



l + (2a-2)-+ 2a J -2a-5 
p 



and 



(44) 



(45) 



(46) 



p — 1 



P p2 ^3 5 



which agree with Eqs. (^TJ) and (|43|). 



(47) 



2.5.2 Inflation near a maximum 

The potential sufficiently near a maximum can be written 

1 



V(<j)) = V [1- 



Therefore, for small 



(48) 



ei = , 6i 



4 
3' 



1 + ^2-1) , 52 = 51, S 3 = 5 ? { 



(49) 
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and the power spectrum and spectral index are 



Vn 



I2ir 2 cj> 2 6 2 
(2tt) 2 2 



r(l) 



7T 



1 + 2a5i + ( 2a z - 4 + — ) 5f 



and 



(50) 



(51) 



which again agree with Eqs. (|4lD and 



2.6 The Spectrum and Spectral Index in terms of the Inflaton 
Potential 



Now we express every quantity in terms of the potential V((f>). Using 3H = V + \<j) , its 



derivatives and H 
Wi = 



2 /2, and defining 



we have 



H 2 

2 



r 
¥ 
i 

i 



Vi 



I! 



V 2 = 2 and V 3 



6 



12 1 9 



i + ~Z4 - ^ - ^ + Zz/ lVl - ^v 2 - ^v 2 

i W 2 + l^ iVi + 4 3 _ 5^2 Vi + 5 ^ 2 + 1 
3 3 9 o 18 9 



2^ 

\ Ul - Vi - jjw? + |w x Vi - ^v 2 - ^v 2 



23, 



17. 



+ - 4^Vi + yWxV^ - -VI + -U X V 2 - - ViV 2 - -V 3 



(52) 

(53) 
(54) 
(55) 

(56) 



U{ - -WiVi + v x 2 + v 2 



^v 2 + \vt - \u,v 2 



2ViV 2 + -V 3 



7 45 

-W 3 - — U\V X + 9WxV 2 - V? + 4W X V 2 - 4ViV 2 - V 3 



Substituting into Eq. fl4"T| ) gives 



Y2-k 2 U x 

'3 2 1 
-a H — a — 
2 3 

: I 2a 2 + -a- 



V 

1 — — 


)Wi + 


77 


ll7T 2 \ 


y + 


8 ) 


37 

Y + 





-2a--) Vi 



W 2 + 



-2a' 



-a + 
3 9 



142 llvr 5 



6 



—a a 

3 9 



TP 



12 



(57) 
(58) 



WiVi 



(59) 



The numerical values of the coefficients in this equation are 2.022245, —2.125941, 1.779141, 
-3.867662, 2.374857, -0.418550, respectively. 
Using 



~ - 1 + -Wi + -Vi + — Wf - -ZYiVi + -Vf + -V 2 



(60) 



ai?=fe 



2 



cilnfc V 6 3 36 1 9 9 9 / 

or substituting into Eq. (fl"3"l) gives 

n^ c (fc) = 1 - 3Wi + 2Vi + (to - W 2 + (-8a - l)W x Vi + jjv 2 + M2a + V : 

/ _ o 13 1867 llvr 2 \ o / „ , 1 595 _ 2 V, 2 ,, 
+ (-12a 2 + y a - _ + _ j U\ + ( 2 4a 2 + -a + — - ll. 2 j ^Vt 

+ (-8a 2 - 6a - ^ + ^) ^V 2 + \v\ + (-6a 2 - 2a - f + 2tt 2 ) U,V 2 
( , 8 28 13tt 2 \ / 2 2 2 vr 2 \ 

The numerical values of the coefficients in this equation are 3.544490, —6.837097, 2.125941, 
-0.804970, 3.621120, -3.800854, -1.247621, 1.119331, 0.418550, respectively. 
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